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Abstract 

How tropical tree species respond to disturbance is a central issue of forest ecology, conservation and resource 
management. We define a hierarchical model to investigate how functional traits measured in control plots relate to the 
population change rate and to demographic rates for recruitment and mortality after disturbance by logging operations. 
Population change and demographic rates were quantified on a 12-year period after disturbance and related to seven 
functional traits measured in control plots. The model was calibrated using a Bayesian Network approach on 53 species 
surveyed in permanent forest plots (37.5 ha) at Paracou in French Guiana. The network analysis allowed us to highlight both 
direct and indirect relationships among predictive variables. Overall, 89% of interspecific variability in the population 
change rate after disturbance were explained by the two demographic rates, the recruitment rate being the most 
explicative variable. Three direct drivers explained 45% of the variability in recruitment rates, including leaf phosphorus 
concentration, with a positive effect, and seed size and wood density with negative effects. Mortality rates were explained 
by interspecific variability in maximum diameter only (25%). Wood density, leaf nitrogen concentration, maximum diameter 
and seed size were not explained by variables in the analysis and thus appear as independent drivers of post-disturbance 
demography. Relationships between functional traits and demographic parameters were consistent with results found in 
undisturbed forests. Functional traits measured in control conditions can thus help predict the fate of tropical tree species 
after disturbance. Indirect relationships also suggest how different processes interact to mediate species demographic 
response. 
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Introduction 

Functional ecology assumes that biological functions scale up 
from the organismic level to higher organization levels [1—4]. At 
the individual level, biological functions can be investigated by 
means of functional traits which are morphological, physiological 
or phenological characteristics that influence individual fitness 
through growth, fecundity and survival [4]. In turn, the 
performance of all individuals in a population integrates into the 
population change rate as assessed in demographic studies, that is 
the outcome of the demographic processes of recruitment and 
mortality. Hence the relationships between the population change 
rate and functional traits are not direct but indirect through 
demographic parameters. Dependencies among population 
change, demographic processes and functional traits are widely 
acknowledged [3,5,6], but have rarely been integrated into a 
consistent model. Here, we evaluate an ecological model [7] 
connecting these variables in a three-level hierarchy (Figure 1), in 
order to identify key traits indicative of species response to 
disturbance in tropical forests. Although the model is general, key 
issues concern the nature of relevant traits and how they relate to 
demographic parameters. 



Among key functional traits that capture most of the interspe- 
cific variation in plant strategies [8], and might thus drive species 
performance [3] , Specific Leaf Area (SLA) stands as a proxy of 
mass-based maximum net photosynthesis, and ultimately relative 
growth rate [9], species size at maturity relates to investment in 
perennial structures and access to light [10], and seed mass 
differentiates species with respect to dispersal and survival in early 
life-stages [1 1]. These three key traits describe the LHS scheme of 
[9] and are indicative of species ecological strategies. Beyond 
these, we consider four other functional traits that relate to critical 
components of species response to disturbance. Wood density 
(WD) correlates with growth potential [12], resource allocation, 
and resistance to damage and embolism [5,13-15]. Leaf nitrogen 
concentration (LNC) reflects the concentration of proteins 
involved in photosynthesis, and thus relates to net photosynthesis 
rate [16]. Leaf phosphorus concentration (LPC) also relates to 
photosynthetic capacity and ultimately growth [17-19]. Lastly, 
leaf carbon isotope discrimination (8 13 C) relates to water use 
efficiency [20], which may be critical in large forest gaps where sun 
exposure can cause drought. 

In tropical forests, relationships between demographic processes 
and functional traits have been mostiy investigated in plots in 



PLOS ONE | www.plosone.org 



1 



September 2014 | Volume 9 | Issue 9 | e1 05022 



Functional Traits and Post-Disturbance Demography 



Performance 
currency 



Demographic 
parameters 



b) 



Population 
Change 
Rate 

X 



Recruitment 
rate 



Mortality 
rate 



Species functional traits 



Specific Leaf Area (SLA) 
Leaf Nitrogen Content (LNC) 
Leaf Phosphorus Content (LPC) 
Carbon isotopic composition (8 ,3 C) 



H 



Maximal diameter (D 
Wood Density (WD) 

S -T Seed Volume (SV) 



Figure 1. Ecological and data models, a) Hierarchical theoretical relationships between Performance currency [3] measured at population level 
and indicative of species performance, Demographic parameters measured at population level and indicative of population dynamics, and Species 
functional traits measured at individual level and indicative of species strategies [4]. b) Data model: variables used to quantify the components of the 
ecological model. Performance is evaluated by the population change rate, k, population dynamics is evaluated by recruitment and mortality rates; 
seven functional traits are considered in relation with the LHS scheme [9]. Gray and white boxes indicate variables measured in disturbed and control 
conditions, respectively. 
doi:1 0.1 371 /journal.pone.01 05022.g001 



natural conditions [12] even if there are some recent attempts to 
link functional traits to disturbance tree responses using long-term 
plots from secondary forests [21,22]. Species mortality rate usually 
decreases with increasing wood density [5,14], adult stature [5], 
and seed size [5]. By contrast, evidence of correlates of recruitment 
rate is scarce in wild tree populations. Large seeds usually favor 
survival in early life-stages, but high recruitment rates are often 
associated with the pioneer syndrome of species producing small 
and numerous wind- or animal-dispersed seeds [23,24], whereas 
recruitment efficiency, that is the recruitment rate adjusted for 
species abundance, decreases with increasing species maximum 
height [25,26] or maximum diameter [27]. All together, these 
relationships demonstrate that functional traits can explain 
demographic parameters in undisturbed conditions to some 
extent, as presented in the model in Fig. 1 . 

In this paper, we hypothesize that traits measured in control 
conditions can explain, directiy or indirectly the demographic 
response to disturbance. To test this hypothesis, we investigated 
the relationships between species traits measured in control 
conditions, and demographic parameters after disturbance due 
to logging operations. When investigating relationships between a 
given trait and a given demographical process, correlation may 
reflect a direct link due to a causative relationship, or an indirect 
link because the two variables are under the influence of a third 
one. Hence, some traits may have direct effects on species 
demography in response to disturbance, while others may have an 
indirect influence. Such trait-based approach allows one to better 
understand the underlying mechanisms of why some are successful 
after disturbance, and others not. It may also initiate a step 
forward into predicting the future of tree species communities 
undergoing disturbance over large areas more easily than by 
settling permanent sample plots. We specifically address the 
following questions: (i) how are population changes in response to 
disturbance partitioned with regard to recruitment and mortality ? 
(ii) how are demographic process rates controlled by functional 
traits ? and (Hi) which traits directly VS indirectly influence the 
demographic rates ? 

The purpose of this paper is to address these questions 
simultaneously in a consistent ecological model. Regression and 
variable selection methods can help address the two first questions. 
The third question however supposes a hierarchical structure 



among variables, with traits having either direct or indirect 
influence on demographic parameters. We used the Bayesian 
Network (BN) framework to investigate the hierarchical structure 
of the ecological model presented in Fig. 1. Other approaches, 
such as the Structural Equation Modelling framework also allow to 
quantify the likelihood of a network of relationships given data. 
However, a key aspect of the BN approach is that the structure of 
the Direct Acyclic Graph (DAG) that represents the network is not 
given a priori but can be inferred from the data. We applied the 
BN methodology to investigate the demographic response of fifty- 
three tropical tree species to disturbance. We used long term 
survey data of population dynamics in logged-over permanent 
sample plots located in French Guiana [27] to quantify species 
population change rates as well as demographic parameters. These 
data were completed by species-specific attributes collected for the 
seven functional traits. After investigating correlative patterns, we 
present the BN framework and apply it to the questions raised 
above. We then discuss the inferred relationships and whether they 
are similar to those observed in undisturbed forests. 

Materials and Methods 

Study site and period 

The study was conducted at the Paracou Experimental site 
(5°18'N, 52°53'W) in a lowland terra firme rainforest of French 
Guiana dominated by the Caesalpiniaceae , Lecythidaceae , Chry- 
sobalanaceae and Sapotaceae families. The site is located on 
private land and rented by the research institute (CIRAD). 
CIRAD is given full permission by the owner to conduct studies on 
the land, including the one presented in this manuscript. No 
specific permission was required for the activities conducted 
during the study (in situ measurements and collections). Field 
collections did not involve any endangered or protected species. 

The climate of the area is equatorial with a dry season from 
mid-August to mid-November and a rainy season, often 
interrupted by a dry period in March. The mean annual rainfall 
and temperature are 3040 mm and 26°C. Soils are mostly shallow 
ferralitic and developed on schists and sandstones that form small 
elliptic hills. 

The experimental setting of the site consists in fifteen 9-ha 
permanent sample plots (PSP), 9 of which were submitted to three 
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logging treatments of increasing intensity (Tl, T2, T3), while 6 
remained control plots in natural conditions (TO). Each 9-ha plot 
contains a core zone of 6.25 ha surrounded by a 25 m wide buffer 
zone. In 1984, all trees with diameter at breast height (DBH) > 
10 cm were localized in the 250 x 250 m core zone of each plot 
and botanically identified. Since then, all plots have been surveyed 
annually to evaluate individual growth, as well as recruitment and 
mortality. Logging operations were conducted during two years 
between 1986 and 1988 at the scale of the 9-ha plots. The 
treatments involved selective harvesting for timber (Tl, T2, T3: 10 
commercial trees. ha -1 DBH >50 or 60 cm), and fuelwood (T3: 
30 non-commercial trees.ha -1 with DBH between 40 and 50 cm), 
and additional thinning of non-commercial trees by poison- 
girdling (T2: 30 trees.ha- 1 DBH >40 cm, T3: 20 trees.ha- 1 
DBH >50 cm; [28]). 

In this study, we focus on the consequences of disturbance on 
population dynamics in terms of indirect effects. Direct effects of 
selective logging and poison-girdling include tree damage and 
destabilization by neighbours loss. These effects were materialized 
by the death of wounded trees up to four years after operations 
stopped. Here, we are interested in the response of species to the 
renewed availability of resources (space, light and soil) after 
disturbance. In order to eliminate the direct effects of the 
treatments, we thus calculated the population change rates and 
demographic parameters starting four years after the end of 
logging operations, over the period 1992-2003. We considered all 
trees > 10 cm DBH occurring in the six most heavily disturbed 
plots, i.e. in treatments T2 and T3. 

Focal Species 

Fifty-three species were selected for this study; they were well- 
documented in floras, unambiguously identified by field workers, 
and abundant enough to supply relevant data (see Appendix SI). 
They represented 28 families with a dominance of Caesalpiniaceae 
(13 species), Sapotaceae (4), Euphorbiaceae (4), Clusiaceae (4), 
Lecythidaceae (3) and Myristicaceae (3), which reflects the local 
composition of the forest. The studied species cover a large range 
of ecological strategies from small and medium-sized heavy shade- 
tolerant species to small light-demanding species and large light- 
demanding canopy species, including three pioneer species. In 
order to characterize these strategies, we estimated species growth 
potential in control conditions. We measured growth potential 
(RDIR for relative diameter increase rate) as the 95'* percentile of 
the distribution of change rate in relative diameter increment 

A In (A) 
At 



D t being the DBH at year t, and At taken here between 2000 
and 2003 [27]. Note that, because there is no direct link between 
population change rates and species RDIRs, we chose to not 
include the latter in the presented bayesian network analysis. 
However, we reported their values in table S 1 in Appendix S 1 and 
performed the network analysis including RDIR for comparison 
purposes (Figure SI). 

Demographic Parameters and Population Change Rate 

We calculated the population change rate over the study period, 
A, as 



where N t(j and N, { are the number of living trees at the beginning 
(1992) and at the end (2003) of the study period of length At. In 
order to take advantage of the yearly surveys, we calculated 
recruitment, defined as ingrowth over 10 cm DBH threshold 
(fREc), and mortality (rjyiOR) rates following [29] as 



Trec, Tmor : 



t' 1 ' 1 n- 



ln(N h -N t() ) 



where ,+ 1 is the number of trees respectively recruited (i.e. 
reaching 1 0 cm DBH) or dying between years / and i + 1 , and N, 
is the number of living trees in year i. We preferred this approach 
to estimate demographic parameters because estimations based on 
the raw comparison of an initial and a final survey [30] cannot 
account for trees both recruiting and dying during the study 
period. We calculated mortality rates excluding trees that had 
been wounded or poisoned during silvicultural operations in order 
to avoid direct effects of disturbance and capture intrinsic specific 
responses to disturbance. Note that although populations change 
in number as a result of the demographic processes of recruitment 
and mortality, the definitions of the corresponding rates do not 
lead to X being a simple combination of the demographic 
parameters. 

Functional Traits 

Leaf traits were measured on leaves collected in three 
undisturbed control plots of Paracou following standardized 
protocols [31]. As far as possible, local measurements in control 
plots have been preferred to values from global trait databases in 
order to obtain locally accurate estimates of the interspecific trait 
variation [32]. Leaves were taken in fully sunlit positions, except 
for three species of small size (Iryanthera hostmannii, Oxandra 
asbeckii and Sandwithia guianensis). Ten leaves in total were 
collected on 3 different trees of each species. Analyses were 
performed on complete leaves, with petiole and rachis, except 
when rachises were too hard to be ground (Carapa procera, 
Jacaranda copaia, Schefflera decaphylla, Simaba cedron). 

We measured the specific leaf area (SLA) as the one-sided area 
of fresh leaves divided by the oven-dried mass. Leaf nitrogen and 
phosphorus concentrations (LNC and LPC, the total amount of N 
and P per unit of leaf dry mass) were measured respectively by gas 
chromatography performed on a Carlo Erba CHN-OS elemental 
analyser, and by the molybdate-blue method, using ascorbic acid 
as a reductant. Leaf carbon isotope composition (8 13 C) was 
measured using a Finnigan DELTA S isotopic ratio-mass 
spectrometer. 

As an estimate of species potential size, we considered the six 
control plots of the site in 2003 for the reference distribution of 
diameters and calculated the maximum diameter (D m ) as the 95* 
percentile of the left-truncated distribution, that is all diameters 
above the threshold of 10% of the species absolute maximum, as 
recommended in [33]. We estimated seed size as the ellipsoidal 
volume calculated from measured dimensions: 



S V = - L x w x cl 
6 



L,w and d being respectively the length, the width and the 
depth of a seed. Seed mass was highly correlated to seed volume 
for a subset of species, but we preferred the latter trait because it 
was available for all studied species. Finally, wood density (WD), 
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the ratio between mass and volume of a piece of wood measured at 
12% moisture content was compiled from a local database and 
literature. Trait data for the focal species are presented in Table 
S2 in Appendix SI. 

Bayesian Network Analysis 

In graph theory, relationships within a set of studied variables X 
can be represented by a Directed Acyclic Graph (DAG). A DAG 
consists in a set of nodes (A,B,C, . . .) and directed edges between 
them representing the variables and the dependencies among 
them respectively: A^B^C for instance. In the Bayesian 
Network (BN) framework [34,35], a statistical model of these 
dependencies is entirely defined by the DAG and an associated 
probability that measures the likelihood of the represented 
dependencies given the data (X). 

In a DAG, the parents of a node n are the nodes that connect 
directly onto n, noted pa(n). When conditioned on its parents, a 
particular node is by definition independent of all other variables 
that are not among the parents. Within each parents-child cluster, 
the dependency of the child node on its parents is thus described 
by a conditional probability distribution, ? (x„\x pa („)) , where x„ is 
the variable associated with node n. The set of distributions for all 
clusters in the DAG allows to define the likelihood of the DAG 
given the data 

p(X|a> n p( Xn \ Xpa{n) ) 



Specifically, the conditional distributions are modelled as 
Gaussian linear models [36]: 

Xn\x plli „)^Af(a + fixx pai „ ) ,a 2 ) (1) 

where a and ji are regression coefficients and a 1 is the conditional 
variance. We thus log-transformed some variables in the analysis 
(X, r REC, r MORand SV) to ensure that their distributions 
approached normality (Table S2 & S3 in Appendix SI). 

Conditionnal dependence or independence among the studied 
variables can be direcdy read from a DAG. In the simple example 
^4<— 5— >C, both nodes A and C depend on B, so that A and C are 
dependent when conditioned on B, i.e. when one considers B as 
fixed [36]. In the example A— >B—>C, node B depends on A, and 
C depends on both A (indirectly) and B (directly), but A and C are 
independent conditionally on B, that is A and C can vary 
independently when B is fixed. Finally, in the example A^>B<— C, 
node B depends on both A and C, but A and C are independent 
because no edge exists between them. In the present study, the 
DAG has an inherent hierarchical structure given by the 
considered ecological model (Figure 1). In order to impose the 
sense of some of the arrows, relationships were specifically banned 
from the network [36] : we did not allow arrows directed from X to 
demographic parameters or traits, neither from demographic 
parameters to traits, or from traits to X. 

We used a score-based approach to infer the structure of the 
DAG from the observed data: the algorithm goes through the 
space of possible networks, scores candidates based on the 
likelihood and finds the best scoring network [36]. We used a 
weakly informative prior for the joint probability, thus assuming 
weak a priori knowledge about the network structure. After a likely 
network was found, we performed an intensive heuristic search, 
which consists in perturbing the DAG a large number of times and 



performing multiple random restarts using the new network. This 
approach insures that a global maximum score has been reached. 

We then performed complementary tests to ensure that the 
inferred multivariate relationships implied by the network were 
statistically valid. We first tested the overall structure of the 
network using Fisher's C test as described by [37]. We then tested 
all inferred relationships implied by the DAG structure ([37,38], 
see Table S4 in Appendix SI) using partial correlation coefficients 

-cab 
rA,B\x = , = 
s/cax x C B X 

where ya,b\x measures the correlation between variables (A,B) 
conditionally on variables in X, and c is the correlation between 
corresponding variables. We used Pearson's correlation coefficient 
to ensure consistency with the linearity hypothesis of the BN 
approach, as well as Spearman's correlation coefficient. All 
inferred conditional independencies were eventually verified. 
Inferred conditional dependencies present in the calibrated 
network that were not significant with regards to both coefficients 
(p>0.05) were removed from the analysis to ensure a conservative 
approach. P-values were corrected for multiple testing. All analyses 
were performed with R [39] using the deal [36] and ggm packages 
for BN and DAG analysis. 

Results 

At the beginning of the study period, population size ranged 
from 19 individuals, in Balizia pedicellaris (Fabaceae) and 
Schefflera decaphylla (Araliaceae) to 1330 in Lecythis persistens 
(Lecythidaceae) for an average population size of 185.2 individuals 
(s.d. = 235.6, see Table SI in Appendix SI). Population change 
rate over the study period varied from -0.01 (Pogonophora 
schomburgkiana Euphorbiaceae) to 0.25 (Schefflera decaphylla), 
while recruitment and mortality rates respectively varied from 
0.4% (Tapura capitulifera Dichapetalaceae) to 15.1% (Schefflera 
decaphylla) and from 0 (Platonia insignis Clusiaceae) to 2.4% 
(Pogonophora schomburgkiana, Table SI in Appendix SI). The 
BN analysis evidenced direct and indirect relationships among 
variables in a network supporting the ecological model in Figure 1 . 
Population change rate, X, direcdy depended on the demographic 
parameters (frec^mor) by construction. Overall, 88% of its 
variance were explained by those parameters, Trec being the most 
influential variable (see arrow width in Figure 2) compared to 

r MOR- 

The recruitment rate Trec was directly related to three 
functional traits that together explained 45% of its variance 
(Figure 2): LPC had a positive effect, whereas wood density (WD), 
and seed volume (SV) had negative effects on Trec- When species 
growth potential (RDIR) was included in the analysis, it appeared 
in a direct relationship with Trec an d the proportion of explained 
variance raised to 54% (Figure SI). The mortality rate was 
explained by a single functional trait, species maximum diameter 
(Dmax, 25%). Among functional traits, LPC was explained by WD 
and LNC (59%). SLA was in turn informed by LNC and D max 
(51%). Lastly, carbon isotope discrimination, 8 13 C, was informed 
by LPC and D max (33%). The relationships between functional 
traits and demographic parameters were only slightly modified 
when we included species growth potential in the analysis. The 
only difference in the structure of the DAG was that the 
relationship between WD and Trec became indirect and mediated 
by species growth (Figure SI): in this model, the recruitment rate 
was then independent of wood density when conditioned on, i.e. 
for given, growth potential and leaf phosphorus concentration. 
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Figure 2. Results of the Bayesian Network analysis. Directed Acyclic Graph (DAG) showing the network of relationships between the 
population growth rate (/), demographic parameters (rREC. recruitment rate, Tmor, mortality rate) and functional traits (see main text for labels). Grey 
(white) indicates variables measured after disturbance (in control conditions). Dashed boxes indicate variables having no parent (explicative variable) 
in the network: SV, WD, LNC and D max . The figure represents a statistical summary of the inferred relationships. Numbers indicate the posterior 
estimates of the model parameters: the intercept a. and the proportion of variance explained by the parents (r; see relation 1) are given in boxes, the 
regression coefficients (/?) are shown on arrows. Arrow width is proportional to the corresponding standardized regression coefficient. Italic numbers 
are estimates of the residual variance (<r 2 ; see relation 1, to be compared with raw variances in Table S3 in Appendix SI). 
doi:1 0.1 371 /journal.pone.01 05022.g002 



The network analysis highlighted four variables without parents 
(exogenous variables; Figure 2): SV, WD, LNC and D max . SV 
appeared on a single path to Trec, while WD, LNC and D mllx 
appeared on forks, each being directly related to two descending 
variables. The positions of these exogenous variables, and 
assuming that errors on them are independent, indicate that these 
variables are independent pairwise (as shown in Fig. 3). Moreover 
their positions as exogeneous variables highlight their influence on 
other relationships. For instance, the position of D max as unique 
parent for Tmor indicates that the mortality rate was independent 
from the other exogeneous variables when conditioned on 
maximum size (Table S4 in Appendix SI). Hence, the correlation 
observed between Tmor and SV (Figure 3) was mediated by D max : 
the two variables appeared independent when conditioned on 
maximum size (Table S4 in Appendix SI): for a given maximum 
size, the mortality rate is independent of seed volume. 

The DAG allows us to highlight independences involving not 
only exogenous, but also endogenous variables (see Table S4 in 
Appendix SI for a list of relationships). For instance, pairwise 
correlation shows that 8 13 C correlates negatively to the mortality 
rate (Fig. 3). However, according to the DAG, the two variables 
are not in a direct relationship and they are in fact independent 
when conditioned on their parents D max and LPC (Table S4 in 
Appendix SI). This finding can be interpreted as mortality rate 
being independent from carbon isotope composition in species 
having similar maximum diameter and leaf phosphorus concen- 
tration. Regarding functional traits, the same reasoning applies to 
the covariation between LPC and SLA, that are independent 
when conditioned on WD, LNC and D max , and SLA and 5 13 C, 
independent when conditioned on LNC, D max and LPC (Table S4 
in Appendix SI). We note that these findings are contextual, they 
derive from the statistical properties of the studied dataset and may 
not be straightforward to interpret biologically. 



Discussion 

Our analysis demonstrates that functional traits measured in 
control plots can explain species demography assessed after 
disturbance. We highlighted both direct and indirect relationships 
between population change rate, demographic parameters and 
functional traits of tropical trees by taking advantage of the 
statistical dependencies among the variables to find the most 
explicative relationships of the hierarchy proposed in a simple 
ecological model (Figure 1). 

Population Change Rate and Demographic Parameters 

After disturbance, recruitment critically controlled species 
demography, while mortality had revealed lower influence 
(Figure 2). Disturbance by logging induced instantaneous demo- 
graphic imbalance in tree populations, because they were either 
directly eliminated by selective logging or poison-girdling, or 
indirectly killed by falling trees or following disturbance-induced 
stresses, implying that most of the final tree populations however 
recovered from reduced numbers [28]. This has resulted in 
increased recruitment [27], as well as increased growth after 
disturbance [40] when compared to populations in control plots. 
Regarding mortality, logging effects were more contrasted when 
compared to control conditions, with decreasing mortality rates 
observed in heliophilous species and increasing rates observed in 
shade-tolerant species, which indicated a negative effect of open 
canopy conditions on survival in the latter (see Table SI in 
Appendix SI). These findings are consistent with previously 
documented disturbance effects on dynamics in logged-over stands 
[41]: dipterocarp forests showed higher recruitment rates in 
treated plots compared to control plots several years after 
mortality rates had returned to pre-disturbance levels. In forests 
of Venezuela, process-based simulations showed that ingrowth 
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Figure 3. Correlation analysis. Relationships between population change rate W) and demographic parameters (r REC , r M oR> after disturbance 
and functional traits measured in control conditions (see main text for labels definition). Panels show respectively variables distribution in the 
diagonal, pairwise connections with lines showing linear regressions in the top triangle, and Spearman correlation coefficients on rank with tests 
results in the lower triangle. Large numbers indicate significant correlations with p-values as follows: * * * <0.001, * * <0.01, * <0.05,~ <0.1. I, 
trec, fMOR, and SV were log-transformed. 
doi:10.1371/journal.pone.0105022.g003 



dominated over mortality after disturbance, leading to positive 
demographic imbalance [42]. We note that the time lag between 
logging operations and the study period, the length of the study 
period and the size threshold used to sample populations all 
determine how populations were sampled and may have 
consequences on the patterns highlighted here. Recruits come 
from germinated seeds or were present as advanced regeneration 
that has grown up to the threshold before the end of the study 
period. Because of their very low inherent growth rate, it is likely 
that recruitment in the more shade-tolerant species should come 
from advanced regeneration. A longer period of censuses would 
allow us to test the hypotheses of (i) a delayed and higher 
recruitment rates for shade-tolerant species together with (ii) a 
concomitant decrease in recruitment of light-demanding, and 
usually small-seeded species. 

Drivers of Demographic Parameters 

About half of the interspecific variability in the recruitment rate 
was best explained by three functional traits, including leaf 
phosphorus concentration, wood density and seed size (Figure 2). 
In disturbed stands of the site, we know that a positive relationship 
exists between growth and recruitment rates [27]. Thus the 
revealed relationships between functional traits such as leaf 



phosphorus concentration and recruitment rate could be mediated 
by species growth potential. Phosphorus is a major component of 
plant growth [43], as it is present in nucleotides implied in 
metabolic pathways (e.g. adenosine triphosphate, ATP). Under 
well-developed vegetation and on ageing soils, weathering and 
biological mobilization decrease phosphorus availability, which 
may in turn limit plant growth and development [44,45] . Nitrogen 
is also known to limit tree growth [43], however N:P ratios 
observed here revealed high values, from 16.9 to 38.5 (26.2 + 5.4), 
which may indicate limitation by P, combined limitation by both P 
and N or even N limitation [45]. Recent evidence however 
confirmed that phosphorus rather than nitrogen controled the 
productivity of Amazonian forests [46]. Here, our network analysis 
evidenced the central role of phosphorus in the system, suggesting 
that species best at mobilizing available phosphorus also better 
recruited in response to disturbance. Unexplained variability in 
recruitment rates proceeds from unobserved causes that are 
sources of variation in juveniles survival and growth, it relates to 
stochasticity in the effects of abiotic and biotic factors on these 
processes. 

Wood density negatively influenced recruitment rate, which is 
consistent with the well-known negative correlation between wood 
density and growth [15,33,47-49]. High wood density implies high 
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stem construction costs, mechanical stability and resistance to 
cavitation [14,50,51], as well as low photosynthetic rates [52]. By 
contrast, low wood density is associated with efficient resource 
acquisition, which allows species to reform population numbers 
rapidly after disturbance. At individual scale, wood density is 
negatively linked to the tree ability to accelerate its growth during 
early life stages [53], that is light-wooded species have the best 
potential to accelerate their growth at intermediate sizes in suitable 
conditions, resulting in a hump-shaped growth trajectory. On the 
contrary, heavy-wooded species with constrained construction 
potential generally show slow [12], but relatively constant, growth 
patterns [53]. 

Seed volume negatively influenced recruitment directiy, sug- 
gesting that regeneration processes were involved in species 
response to disturbance: high disturbance intensity tended to 
preferentially favor small-seeded species. However, as found 
elsewhere, we note that correlations between demographic 
parameters and seed volume were weak [5]. In undisturbed 
tropical forests, seed size was negatively correlated to mortality [5] . 
Although not causally linked to mortality of large trees, this 
correlation can be interpreted as a consequence of the adaptive 
trade-off between tolerance to low resource supply and growth 
[1 1,23]. Large seeds often occur in species with attributes related 
to shade-tolerance, whereas pioneer and light-demanding species 
disperse overall smaller and more numerous seeds [54]. 

Previous studies found that mortality rates were negatively 
correlated with wood density in undisturbed tropical forests 
[5,33,55,56]. In our case, wood density was unrelated to mortality 
rate measured after disturbance. We suspect that logging-induced 
disturbance has blurred connections between the mortality rate 
and wood density. During the study period, light-demanding 
species, which are also light-wooded, survived better in disturbed 
than in control conditions (Figure 2) due to higher light availability 
[27], while disturbance had lower effects on the mortality rate of 
more shade-tolerant species. It also appears that correlations 
between wood density and demographic parameters might be 
stronger when considered over larger spatial scales (but see [6,56]) 
leading to wider variability in wood density in relation with, for 
instance, variation in soil fertility [57,58]. 

Species of high stature showed significantly lower mortality rates 
than smaller ones, consistently with results observed in undis- 
turbed forests [27,56]. We suspect that disturbance reinforced this 
relationship in logged-over plots: species of small stature are often 
shade-tolerant and, although some of them benefited from canopy 
opening [40], the majority probably suffered from sun exposure. 
In logged forests, gap creation increases lateral lighting and 
temperature along the vegetation profile. Such conditions have 
also likely impacted species differentially according to their ability 
to optimize water use. Interestingly, mortality also varied with 
enrichment in heavy carbon isotope, 8 13 C. In sunlit leaves, 5 13 C 
relates positively to water use efficiency (WUE; [59]), which 
measures the ratio between biomass gain through CO2 assimila- 
tion and water loss through evaporation. At Paracou forest site, 
water availability is the main climatic driver of tree growth [60]. It 
has been shown that tall tree species preferentially extract water 
from layers below 100 cm depth while shorter tree species show 
broader variations in the depth of water uptake [61]. Here, 
mortality decreased in species with increasing values of carbon 
isotope composition (higher WUE), but the correlation was 
indirect: the relationship appeared mediated by species stature 
and leaf phosphorus concentration. For a given maximum size, 
and for species with similar LPC, susceptibility to drought as 
indicated by carbon isotope composition, was not found to 
influence tree mortality. 



Our network analysis evidenced the role of four functional traits 
as independent sources of interspecific variability. More specifi- 
cally, leaf nitrogen concentration, wood density, maximum 
diameter and seed volume appeared as exogeneous variables, that 
is, variables not explained by other variables considered in the 
analysis. This finding constitutes an indirect validation of the leaf- 
height-seed (LHS) scheme which considers three major indepen- 
dent axes of variation in plant ecological strategies [9] . Among 
those four exogeneous variables, the proxy for the "leaf axis" of 
the LHS was leaf nitrogen concentration. Its relationships with the 
demographic parameters however were weak for recruitment and 
only marginal for mortality. Furthermore, SLA correlated better 
with LNC than LPC, LNG was more strongly correlated to LPC 
than SLA, and LPC appeared more strongly correlated to 
recruitment than LNC. Altogether, these relationships resulted 
in LNC acting as an independent axis of variation, LPC being a 
direct driver of recruitment and SLA having no direct influence on 
species response to disturbance. Wood density and maximum 
diameter characterized two different aspects of the "height" axis 
[12] in relation with resource allocation and conservation (WD) 
and resource acquisition (D max ). Eventually, seed volume 
independendy positioned species along the seed axis of plant 
ecological strategies. 

Relationships among Functional Traits 

We found close positive associations among three major traits 
involved in the so-called "leaf economics spectrum" (LES; [16]): 
SLA, LNC and LPC. These associations are indicative of 
coordination among the different processes leading to C02 
fixation by leaves, as SLA relates to light capture, LNC to 
carboxylation capacity, and LPC to energy transfer during 
photosynthesis. High values of these three traits lead to high rates 
of photosynthesis and respiration [16,62]. 

The relationships between leaf carbon isotope composition, 
8 13 C and the other traits of the LES are not straightforward to 
interpret. Relationships with leaf nutrients have not been found to 
be consistent across studies [63-65]. As an integrative proxy of 
WUE, 8 13 C is strongly driven by stomatal conductance and 
maximum photosynthetic rate. Here, leaf 8 13 C was found to be 
positively related to LPC (directly), but not to LNC, which may 
indicate a control of WUE by species photosynthetic activity 
rather than water flows. If one assumes low variance in stomatal 
conductance across species, these findings suggests that higher 
LPC implies higher 8 13 C because of lower concentration of C02 
in leaves due to higher photosynthetic rates [66]. Such a 
relationship was not observed with nitrogen, which also suggests 
limitation of photosynthesis by P rather than N. 

A low water use efficiency (low leaf 8 13 C values) appears to be 
often related to a high SLA [63-65,67,68], but the lack of a 
relationship with SLA has already been found elsewhere [69,70]. 
Clearly, a better understanding of how leaf 8 13 C relates to traits of 
the LES is required, which implies a thorough analysis of how 
carbon and water economies are coordinated within the leaf. We 
also found that leaf 8 13 C increased with species maximum 
diameter, a surrogate of maximum height. This pattern has been 
reported elsewhere both within [71,72] and across species [20]. 
Gravity and xylem path length resistance induce a reduction in 
water potential at tree tops with increasing species height. This 
reduction is usually compensated for by decreased stomatal 
conductance leading to C 13 -enrichment and increasing 8 13 C 
[72,73]. 
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Conclusions 

The present study reveals the potential of functional traits 
measured following standard protocols to explain and predict 
species performance in response to disturbance and logging 
operations in particular. Regarding recruitment, 45% of the 
variability in the process rate across species could be explained by 
three traits only, WD, SV and LPC, whereas for mortality one 
trait only, D max , significantly explained the process rate up to 
25%. These results are consistent with other studies considering 
functional traits as drivers of demographic processes (e.g. [6]). The 
discrepancy between the two processes partly results from the 
higher part of stochasticity in mortality in tropical forests. We used 
a novel approach to unravel indirect relationships showing a 
hierarchy of traits being more or less closely related to 
demographic parameters at population level. However, we 
acknowledge that the modelling framework used here did not 
allow interactions among traits to drive the demographic 
parameter values, although interspecific variation in demograph- 
ical rates in nature is shaped by syndromes of functional traits 
[74,75]. Interestingly four traits among those tested were 
independent of all others and therefore form a minimal set of 
predictors of species response to disturbance. LPC was however a 
better predictor of recruitment rates than LNC. These results 
evidence that local databases of functional traits can help 
understand and predict the response of forest communities to 
anthropogenic disturbance (see also [76]). Understanding tropical 
forests response to disturbance, including anthropic, is an 
important challenge to scientists which ideally requires long-term 
data collected in permanent plots with different disturbance 
regimes. These plots are often costly to manage, and cannot 
sample each particular "site x forest type" combination neither 
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